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We describe an application of generalised unitarity to the computation of one-loop amplitudes with massive 
external fermions. We present analytic results for the cut-constructible parts of the leading colour contributions 
to the all-plus helicity configuration of the tigg amplitude. Using a special choice for the helicity basis of the 
massive fermions we are able to obtain extremely compact analytic expressions. In particular we describe how 
one can fix the divergent contributions from tadpole and wave-function renormalisation using universal UV and 
IR behaviour. 



1. Introduction 

The Large Hadron Collider, due to start later 
this year, will require accurate QCD predictions 
for successful data analysis. New TeV scale 
physics is likely to be associated with the pro- 
duction and subsequent decay of heavy particles 
and hence background signatures will come in the 
form of multi-jet final states. As quantitative 
QCD predictions only begin at NLO the main in- 
gredient for the computation of such backgrounds 
are multi-particle one-loop scattering amplitudes. 
Traditional Feynman methods struggle to cope 
with the enormous growth in the number of di- 
agrams with the number of external legs and so 
many processes remain uncomputed. As an al- 
ternative, on-shell methods provide an extremely 
efficient tool by keeping only the physical degrees 
of freedom, simplifying calculations. For a review 
of the subject we refer the reader to reference [1]. 

Unitarity based techniques [213] have been used 
successfully for many years to compute compli- 
cated multi-particle and multi-loop gauge theory 
processes. Recent developments utilising complex 
momenta have lead to new recursive techniques at 
both tree level [4] and one-loop [5\ as well as im- 
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proved multiple cutting techniques for logarith- 
mic terms IBl. 

A general one-loop amplitude can be written 
in terms of a basis of scalar box, triangle, bub- 
ble and tadpole functions together with rational 
terms missed by standard four-dimensional cuts. 
Because the quadruple cut completely freezes the 
loop integration, determination of the box coeffi- 
cients is a completely algebraic procedure. For 
the lower-point integral functions Ossola, Pa- 
padopoulos and Pittau have shown how to sys- 
tematically generate subtraction terms to find al- 
gebraic procedures for the remaining coefficients 
[7]. By using a special complex parametrisation 
of the cut momenta, Forde has demonstrated that 
this method can also be used to generate compact 
analytic expressions without subtraction |8j . Ini- 
tial numerical implementations, such as BlackHat 
[5] and Rocket |10lll| . are making the first steps 
towards a NLO event generator. 

In this paper we demonstrate the power and di- 
versity of these new tools by computing the cut- 
constructable contributions to top pair produc- 
tion via gluon fusion. We derive compact ana- 
lytic expressions in agreement with much length- 
ier expressions obtained from previous Feynman 
diagram based computations [12]. This builds 
upon a growing body of work generalising uni- 
tarity methods to massive processes [13 14 15]. 

Throughout this paper all amplitudes are con- 
sidered to be colour-ordered helicity amplitudes, 
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Figure 1 . Quadruple cut for box coefRcients with 
all momenta out-going and loop momentum flow- 
ing clockwise. 



using the standard spinor-helicity formalism to 
write the amplitudes in terms of spinor products. 
We accommodate massive momenta into the for- 
malism by decomposing them into two massless 
momenta using an arbitrary massless reference 
vector, 77 [16]: 



p = p + -;— J— — ry. 



('7|p|'7] 



(1) 



This then allows the definition of external fermion 
wavefunctions as follows [17I18J . 



u±{q, m) 
v±{q,m) 



{i^ + m)\riT) 
{q'±\f]T) ' 

(^-to)|77t) 
{q'±\vT) ■ 
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2. Integral Coefficients with massive prop- 
agators 

In this section we outline the method for direct 
extraction of integral coefficients with arbitrary 
internal and external masses |8|13j . Throughout 
the next section we take all external momenta to 
be out-going and all loop momenta to circulate in 
a clockwise direction. 

2.1. Box Coefficients 

For the coefficients of the box integrals, the on- 
shell constraints on the loop momenta freeze the 
loop integration resulting in a purely algebraic 



procedure for the evaluation of the quadruple cut 
[6]. Here we use a loop momentum parametrised 
by 



h = aK4^ + bKi^ + c\Ki^)[Ki^\ + d\Ki^)[K4^\, 



where 

b 1K4 - S4K1 
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One can find the values oi a,b,c,d from the con- 
straints, 



S={ll ^ml,{h^Kif =ml, 

{h -Ki- K2f = ml, (h + K^f = ml]- 



explicitly this gives us 

5*1 5*4 + 75*1 
7^ - S'iS'4 



d — — [ ah 

c V 7 



(8) 



S'i54-l-7S'i 
b= ^-^, (9) 



7^ - S'iS'4 



(10) 



and c as a solution to the quadratic equation: 



+ c (a{Ki'\K2\Ki'] + b{K^'\K2\Ki'] - S12) = 0. 

(11) 

Here the mass dependence is determined by 54 = 
5*4 + m,\ — to|, Si — Si -\- m^ — TOj and iS'12 = 
S'2 -I- 2Ki ■ K2 + m2 — m.^. We label the complex 
solutions for the loop momentum as li , corre- 
sponding to the solutions c± of eq. (fTT|) . The 
value of the box coefficient is then given simply 
as the product of four tree amplitudes evaluated 
at the values, li , of the loop momentum: 



^m 1 m 1 7712 ^Ti3 m4 
^4;Ki\K2\K3 2 



'- J2 AiA2A,A,{l'[). (12) 



cr=± 
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Figure 2. Generic triangle cut diagram 



2.2. Triangle Coefficients 

For the triple cuts we are left with a non-trivial 
one-dimensional integral after the on-shell condi- 
tions have been imposed. Following the formal- 
ism of Forde [8] we can use simple complex anal- 
ysis to isolate the scalar triangle coefficients, and 
again reduce the problem to a purely algebraic 
one. We first write down a generic triple cut in- 
tegral as 



loop momenta, 1°' , are the complex solutions to 
the on-shell conditions, 

S = {lj= ml {h-Kif = ml, {h+K:if = ml). 

(16) 

Because the second term in the above equation 
has no t dependence in the numerator and has at 
least one additional propagator it must be asso- 
ciated with the scalar box coefficients which have 
already been determined. Therefore the only con- 
tribution to the scalar triangle is the boundary 
behaviour described by the large t polynomial, 
Inf t . For the given choice of parametrisation, the 
integrals over positive and negative powers of t 
vanish: 



dtJtt'' = VneZ/{0}. 



(17) 



We are then only left with the i" component 
which forms the complete coefficient, 



^rnim2m3 
^i-Ki\K2 



-lY.^nU[A,A,A,{r,) 



(18) 
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d^hl[6{l^-m^^)A,A2A3. (13) 



i=l 



The loop momentum is chosen in a similar way 
to the box case but now has a free a complex 
parameter, t: 

h = aKs^ + bKi^+t\Ks^)[K,^\ + j\K,^)[K3^\, 

(14) 

where Ki , K3 are defined analogously to eqs. 
([5][7]). The triple cut integral can then be written 
in terms of its pole structure in the complex t 
plane, 

^3 ^ 5J/"^*^*E (l^idAiA.Asill)] 
Rest^tMiA2A3il^) 



t-U 



(15) 



2.3. Bubble Coefficients 

The bubble coefficients can be extracted from 
the double cut using a similar analysis to the 
triple cut case, as considered in the previous sec- 
tion. In this case two non-trivial integrals remain, 
which we choose to be parametrised by y and t. 
The basis for the loop momentum in this case re- 
quires the introduction of an arbitrary massless 





where the Inft operation encodes the polynomial 
behaviour on the boundary of the t contour. The 



Figure 3. Bubble cuts with triangle subtraction 
terms 
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vector X for which the final coefficient is indepen- 
dent, 

h=yK,^ + ail~y)x + t\K,^)[x\ 



+ 



where 



-\x)[Ki 
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1 = (xl-^ilxl- 



(19) 



(20) 



The coefficients a and h are solutions to the on- 
shell constraints, 



{ll=ml{h-K^f=ml}. 



(21) 



We find that there are non-vanishing integrals 
over both y and t which complicates the proce- 
dure and requires both pure bubble and triangle 
subtraction terms to be included. The final for- 
mula for the coefficient is relatively simple how- 
ever, 






\v^^Yi,t« 



-tlnU[lniy[AiA2ilT, 

-^ E E Inft[Res,[^iA2A3(/?,y±)]|,.^T., 
{K2}y=y± 

(22) 

where y± are the two solutions to the additional 
triple cut constraint {{li + K^)"^ — ml} . The ex- 
plicit values for these non- vanishing integrals can 
be found in references [8113119] . although some 
care must be taken when dealing with the differ- 
ent momentum conservation conventions. 

3. Application to top pair production 

The four-point amplitudes for top pair pro- 
duction through gluon fusion has were computed 
long ago using Feynman techniques [20j and pre- 
sented analytically more recently [12]. As a test 
of our method for extracting integral coefficients 
for arbitrary internal and external masses we re- 
compute the cut-constructable parts of the colour 
ordered primitive amplitudes, A]^ and A4 , de- 
fined by: 



i(i) 



m, 



A^^^(lQ,2,3,4Q) = Ai"'(lQ,2,3,4Q) 
+ -r^^i (1q: 2, 3, 4q) -t- fermion/scalar loops 



rational terms. 



(23) 



The left-moving primitive amplitudes are those in 
which the fermion turns left around the loop and 
the right-moving amplitudes where it turns right. 
The full integral basis is shown in figure [H In 
terms of the integral basis, the cut-constructible 
parts of the primitive amplitudes are: 



AWn o o A \ ^OOOm rOOOm 

^4 (-Lq,^,3,4qj = G4 74 



OOm 7-OOm 



COniO rOrnO 
3;23 ^3 



'-^3;12''3 +'-"3;23^3 

/-.OOO tOOO I r^Om j-Om , /-(OO rOO 
"-^3:23^3 T" ^2:12^2 "*" "-"2:23^2 



CrCi / ^1 



(24) 



mmmO rmmmO 



4''klQ,2,3,4Q) = C4"™"%" 

I ^mmO 7-mmO i^ /-imniQ TmmO i^ ^mOm jraQm 
+ '-"3;23 -"S 



-'3;23 ^3 
''2:12-'2 



-'3;23 ^3 



L'2:12-'2 "T "-^2:23 -'2 



+ 



ctC[ ( fi 



(25) 



The compact tree level input for each coefficient 
can be generated through on-shell recursion rela- 
tions 4 . The procedure for computing the co- 
efficient is easily automated and has been done 
both numerically and analytically. The analytic 
form of the final coefficient strongly depends on 
the various choices for the reference momenta of 
both internal and external particles. In order find 
the most compact representations it is convenient 
to begin with the most general representation and 
use complex reference momenta. 

Here we present the full results for the 
all-plus helicity configuration, the remaining 
helicity amplitudes will be presented else- 
where. In this simple case the coefficients 

/^OOO r^vamm /^OmO y^OOm y^OOrn /^OO onrl r^rnm 
"-'3;23! '-^3;23 ' "-'3;23 ' '-'3:12 ' L'3;34 : "-'2;23 '^'^'^ '-'2;23 i 

all vanish leaving only 7 out of the 14 non-zero. In 
fact we can further simplify the solution by mak- 
ing a special choice for the reference momenta 
771,774. Choosing 771 = 7?4 = P2 also makes the 
(jooom coefficients vanish and the 6 non zero co- 
efficients are: 



C4"™"°(l+,2+,3+,4+) = 

7r/73(2|i|2][32](2(2|l|2]-h523) 
2(21^)(24'')(3|1|2] 



(26) 
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Figure 4. The scalar integral basis for the ttgg amplitude. The red dot-dashed lines represent massive 
fermions whereas the solid black lines represent gluons. The internal mass labels appear as superscripts, 
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im3(2|l|2][32] 





2(21^)(24^)(3|1|2] 
C2";T2(l5,2+,3+,4+) = 

im3(2|l|3][23] 
(21^)(24^)(2|1|2]2 



C0T2(U,2+,3+,4+) 



(27) 



(28) 



(29) 
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All of these expressions agree numerically with 
expressions extracted from the Feynman calcula- 
tion of ref. fni. 



3.1. Fixing the tadpole and wave-function 
renormalisation coefficients 

The final step in the calculation is to fix the co- 
efficients of the remaining - terms coming from 
tadpole and wavefunction renormalisation contri- 
butions. Although it may be possible to evalu- 
ate these terms using unitarity cuts, a number 
of subtleties arise, as has been discussed in a re- 
cent numerical study [21j . Instead we find it more 
practical to use the universal UV and IR factori- 
sation properties |22|23j . 

Here we present a simple argument which is 
sufficient to find the remaining cut-constructible 
terms for the tigg amplitude. A complete analysis 
will appear elsewhere [19]. The first step is to 
notice that the pure - poles must be proportional 

to the tree level amplitude, A\ . Since we know 
the analytic forms for the bubble integrals, and 
their coefficients this allows us to write down the 
coefficients of the remaining \og{vn?) terms. From 
the basis for the primitive amplitudes, eqs. ([241 - 
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, we see these log(TO^) terms can be completely 
determined from the coefficient of the pure - pole, 



Ol — ayi4 <-^2;12 ^2:23 



C{ = a' A 



(0) 



^2:12 



^2;23- 



(32) 
(33) 



The coefficients a and a' are simply derived from 
the well known UV behaviour of the massless am- 
plitudes and the recent small mass factorisation 
of Moch and Mitov [23] . The factorisation of the 
massive loop amplitude, A^^'(m), is given by, 

n 

A^'\m) ™-^° A(i)(0) + J2 4^'^'^^A^"H0) (34) 



From the known pole structure of the massless 
amplitude, ^(^^(0) , we find a contribution of 
-^H^ [24125] to the 1/e coefficient. It is then 
straightforward to read off an additional factor of 
—Cpng from Zj for each heavy quark appear- 
ing in the amplitude g Putting the two pieces of 
information together gives us: 



a — —a 



2' 



(35) 



which matches the known result [T^ . 

4. Conclusions 

We have shown that using the unitarity for- 
malism of Forde [S] , and its generalisation for ar- 
bitrary masses [13], provides an efficient method 
for calculations of complicated QCD processes. 
We presented a simple example of a helicity am- 
plitude for the process of top pair production 
through gluon fusion while we refer the reader 
to [in] for a complete analysis. 

We foresee that future developments to include 
rational contributions [10I26I9J should permit ap- 
plications to much needed processes of top pro- 
duction in association with many jets in the near 
future [5T]. 



^ There will also be gluon self energy corrections, Zlf , but 
they appear in heavy quark primitive amplitudes not con- 
sidered here. 
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